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Abstract 

Let n > 3 and JJ„ := circ(Ji, J 2 , . . . , J„) and 2„ := circ(j , ji, ■ • ■ ,jn-i) 
be the n x n circulant matrices, associated with the nth Jacobsthal num- 
ber J n and the nth Jacobsthal-Lucas number j„, respectively. The de- 
terminants of J n and 1 n are obtained in terms of the Jacobsthal and 
Jacobsthal-Lucas numbers. These imply that J n and l n are invertible. 
We also derive the inverses of J n and l n . 



1 Introduction 

The n x n circulant matrix C' n := circ(co,Ci, 
numbers cq, . . . , c„_i, is denned as 



, c„_i), assoicated with the 



c ci 
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(1) 



Circulant matrices have a wide range of applications, for examples in signal pro- 
cessing, coding theory, image processing, digital image disposal, self-regress de- 
sign and so on. Numerical solutions of the certain types of elliptic and parabolic 
partial differential equations with periodic boundary conditions often involve 
linear systems associated with circulant matrices [9-11]. 
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The eigenvalues and eigenvectors of C n are well-known [14]: 

71-1 

Aj = ^2c k aj 3k , j = Q,...,n-l, 

k=0 

where w:=exp(^p) and i := and the corresponding eigenvectors 

Vj = . . . , 6>-^") T , j = 0, . . . ,n — 1. 

Thus we have the determinants and inverses of nonsingular circulant matrices 
[1,3,4,14]: 

7i—l n — 1 

det(C„) = n^E^)' 

j=0 fc=0 

and 

Cn 1 = circ(a 0) ai, . . . ,a n _i), 

where a,j-—-^ J2k=o Afe^ - * 1 - 7 ' , and r = 0, 1, . . . , n — 1 [4]. When n is getting large, 
the above formulas are not very handy to use. If there is some structure among 
Co, . . . , Cn—i, we may be able to get more explicit forms of the eigenvalues, deter- 
minants and inverses of C n . Recently, studies on the circulant matrices involving 
interesting number sequences appeared. In [1] the determinants and inverses of 
the circulant matrices A n = circ(Fi, F%, . . . , F n ) and B„ = circ(Li, L%, . . . , L n ) 
are derived, where F n and L n are the nth Fibonacci and Lucas numbers, re- 
spectively. In [2] the r-circulant matrix is defined and its norm is computed. 
The norms of Toeplitz matrices [13] involving Fibonacci and Lucas numbers are 
obtained [5]. Miladinovic and Stanimirovic [6] gave an explicit formula of the 
Moore-Penrose inverse of singular generalized Fibonacci matrix. Lee and et al. 
found the factorizations and eigenvalues of Fibonacci and symmetric Fibonacci 
matrices [7]. 

When n > 2, the Jacobsthal and Jacobsthal-Lucas sequences { J„} and {j n } 
are defined by J n = J n -i +2J„_2 and j n = jn-i + 2jn-2 with initial conditions 
Jo = 0, Ji = 1, jo = 2, and ji = 1, respectively. Let J n := circ(Ji, J2, . . . , J„) 
and Xi '■= circ(jO) ji, ■ ■ ■ ,jn— l)- The aim of this paper is to establish some 
useful formulas for the determinants and inverses of J ra and 1 n using the nice 
properties of the Jacobsthal and Jacobsthal-Lucas numbers. Question: How 
about eigenvalues? Matrix decompositions are derived for J„ and 1 n in order 
to obtain the results. 

2 Determinants of J n and 2 n 

Recall that J„ := circ(Ji, J2, . . . , J n ) and 2 n := circ(jo,ii, . . . , jn-i), i.e., where 
Jfc and jk are the fcth Jacobsthal and Jacobsthal-Lucas numbers, respectively, 
with the recurrence relations Jk = Jk-i + 2Jk-2, jk = jfc-i + and the 

initial conditions Jq = 0, J\ — 1, jo = 2, and j\ = 1 (k > 2). Let a and ft be 
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the roots of x 2 — x — 2 = 0. Using the Binet formulas [8, p. 40] for the sequences 
{ J n } and {jn}, one has 



-[2» -(-!)"] 



and 



Jn = a 

Theorem 1 Let n > 3. Then 



, n i on 



P n = 2 n + (-1) 



(2) 
(3) 



det(J„) - (1 - J„+i)"- 2 (l - J n ) + 2^ [J fe (l - J„ +1 ) fe - 1 (2J„) n - fc - 1 ] . (4) 



k=l 



Proof. Obviously, det(Js) = 20. It satisfies ([tj). For n > 3, we select the 
matrices P n and Q n so that when we multiply J n with P n on the left and Q n on 
the right we obtain a special upper triangular matrix that have nonzero entries 
only on the first two rows, main diagonal and super diagonal: 



and 



P n : = 



1 

-10 

-2 





1-1 
1-1-2 








2 J, 




l-J r 


+ 1 


2 J, 




1-J r 


+ 1 


2 J, 





n— 3 



2J„ 



1 -J„+i 
1 





1 

-1 

-2 









-1 



-1 



(5) 
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Notice that we have the following equivalence: 



J^n^nQn 










1 fn Jn 


— J-n-1 




-Ja 


-Js 


9n Jn 1 


— 2 J n -2 


— 2 J n -3 


-2J 3 


-2J 2 


2t/ n 


Jn+1 — 1 

2J n 


Jn+1 ~ 1 







2 7 



2(/ n Jn-\-l 1 



and S n is upper triangular, where 

n-i 



fc=i 



-fe-i 



ra-2 



.9?i : — i — j n + 2 j w _fc_i f - 



fc=l 



Then we have 



Since 



and 



det(S„) = det(P„) det(J„) det(Q„) = (2 J„)"- 2 . 9 „. 

1 n = 1, or 2 (mod4) 



det(P„) 



2J„ 



det(Q„) = { V 1 -W. „_ 2 

2J„ \ 
.1- J 

for all n > 3, 

det(P„) det(Q„) = 

and ([|) follows. ■ 
Theorem 2 Lei n > 3. TTiera 



1, n = or 3 (mod4), 

2 

, n = 1 or 2 (mod4) 



= or 3 (mod4), 



2*7*1 



i - Ji 



»+i 



n-2 



detp n ) = (2-j n )"- 2 (4-. 7 „_ 1 ) + E [( 2 J'fc - Jfe-i)( 2 - in) fe " 2 (l + 2j„_i) n - fc ] . 

(6) 



fe=2 
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Proof. Since det(3 3 ) = 104, H 3 satisfies (^|). For n > 3, we select the matrices 
K n and M n so that when we multiply 1 n with K n on the left and M n on the 
right wc obtain a special upper triangular matrix that have nonzero entries only 
on the first two rows, main diagonal and super diagonal: 



K n :~- 



1 

2 

-2 














1 
1 -1 








-1 

-2 





1 

1 -1 

-1 -2 



(7) 



and 



M n := 







l + 2jn-l 

2- On 



n— 3 
n— 4 



l+2j„-] 
2-Jn 
1 






-1 

-1 



We have 



U n = K n l n M n 

[2 y ; 



-Jn-1 



-Jn-2 



Un \jn-l — JO \jn-2 — in-1 
1 + 2j„_l J„ - 2 

1 + 2j„_i 



-J3 -J2 

\h - k \h - h 



jn - 2 

1 + ^ jn-1 jn - 2 

1 + 2j„_i 



and U n is upper triangular, where 



1 

2 

n-l 

E 

k=l 



(4-i„-i) + 5^(2i fc -i fc _ 1 ) 



fc=2 



1 + 2j„-i 
2 - j n 



i-fe-1 



1 + 2jn 

2 - J* 
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Then we obtain 

det(t/„) = det(Jf n )detp n )det(M n ) = 2(1 + 2 Jn _ 1 ) n ~ 2 y n . 

Since 

det(i^„) 



1, n = 1 or 2 (mod4) 
— 1, 77, = or 3 (mod4), 



and 

2-j„ 

\ n-2 

l+2j„ 
2-j„ 



det(M„) 
for all n > 3, 



77 = 1 or 2 (mod4) 
, 77 = or 3 (mod4), 



l + 2j»_ 1 N " 
2 - i„ 



det(Jf n ) det(M„) = 
and we have ([|) . ■ 

3 Inverses of J n and H n 

We will use the well-known fact that the inverse of a nonsingular circulant matrix 
is also circulant [14, p.84] [12, p.33], [4, p.90-91]. 

Theorem 3 The matrix J„ = circ(Ji, J2, . . . , J n ) is invertible when 77 > 3. 

Proof. From Theorem g, dct(J 3 ) = 20 ^ and dct(J 4 ) = -400 ^ 0. Then 
J3 and J4 are invertible. Let 77 > 5. The Binct formula for Jacobsthal numbers 
gives J n = - ~^ , where a + /3 = 1, a/3 = —2 and a — (3 = 3. Then we have 



n (, M \^ 7 kr — k ST^ I a P \ , kr-k 1 ST^ ( „r or\ , kr-k 

9(u ) = l^ JrUJ 3 ) =32^( a ~P) UJ 

r—1 r—1 ^ ' r— 1 



(1 - a n ) (3(1 - p n ) 



, (1 - aw fc ,l - (3uj k ^ 0) 



1 — auj k 1 — /3w fe 
1 ((a- 13)- (a n+1 - (3 n+l ) + a(3uj k (a n - f3 n ) 



3 \ 1 — aw' — /3w fe + a(3uj 2k 

1 — Jn+l — 2J n UJ k 



k o, ,2fe ' 



1 - w fe - 2lj 



k = 1,2,. ..,77- 1. 



If there existed u) k (k = 1, 2, . . . , 77 — 1) such that g(to k ) = 0, then we would 
have 1 - J n+ i - 2 J n u k = for 1 - uj k - 2w 2fe ^ 0. Hence w' 1 = 1- 2 j" +1 . It is well 
known that 

(2kiri\ { 2kir\ . /2/c?t 

= ens -I- i sm 



wWp(— j=cosl— j+isinl — ) (8) 



G 



where i := \/— T. Since uj k = 1 2 j" +1 is a real number, sin (^f-) = so that 
ui k = — 1 for < < 2ir. However u — — 1 is not a root of the equation 
1 — J n +i — 2J n u = (n > 5), a contradiction, i.e., g(u! k ) ^ for any u> k , where 
k = 1, 2, . . . , n — 1, n > 5. Thus the proof is completed by [1, Lemma 1.1]. H 

Lemma 4 Lei A = (fly) be the (n — 2) x (n — 2) matrix defined by 

2t/ n , z = j 

J n + l - 1, J = 1 + 1 

0, otherwise. 



Then A 1 = (a^ ) is given by 

( (i-.j n+1 y- 



J > « 



0, otherwise. 



Proof. Let B = = 1 . Clearly by- = ^fc=i a ik<i k j. When i = j, we 
have ^ 

ba 2t/ n * 1. 

If j > i, then 

n-2 



fe=l 



(1- Jn+l)^" 1 , OT (1- Jn+l)^ 



~ (J " +1_1) (2J„)^ +2J " (2J„)i-*+i 
similar for j < i. Thus AA^ 1 = I n -2- M 



Theorem 5 Let the matrix J n 6e JJ„ := circ(Ji, J2, . . . , J n ) (n > 3). Then the 
inverse of the matrix J„ is 



where 



mi 



m 2 



m 3 



J n 1 = circ(m 1 ,ni2, . . . , m n ) 
Jn+i + (1 - 2J„_i)g„ - 1 



JnQn 



(1 - J„ - gr n ) 



(1 - Jn+l) 
(2J„) r 



n— 3 



2^J„_, 



(1 - J„+i)"- fc - 2 



*;=2 



(2Jn) 



i-fe-1 
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_1_ ( {(X -Jn~ 9n)(Jn+2 ~ 1) ~ 4J n J n _ 2 ](l - Jn+l)"" 4 
9n \ (2J„)- 2 



(1- Jn ^ 



'n+1, 



1)1; 



fc=l 

1 (1-J n+1 )" 



(2J„) S 



(1 ' J "^}^" +L ' 4) - 2 + J - a(1 J ; J " +l) ) 

Jn-2(l-J n + l) 



n is even 



/or g n = 1 - J„ + 2 J]fc=i Jn-k-i ( anrf i = 5, 6, . . . , n. 

Proof. Let 



1 ,/n ^ (Jn 1) ~l~ Jn Jn — 1 2 Q ^ t/n- 



1 










J„-l 

1 





2J„- 



n—k—1 



J 3- 2 t J2 

2 J - 2 

5,1 







1 

and g n = 1 - J n 



and G = diag(l,g„) where /„ = £)fc=i Jfc+i (^ 1 _ 2 j^' +1 
(i^fcr) ■ Thcn we can writc 

Pnj nQnRn — G ffi .A 



where G© A is the direct sum of the matrices G and A. Let T n — Q n R n . Then 
we have 

JJ™ 1 = T„(G _1 © A~ 1 )P n . 

Since the matrix J„ is circulant, its inverse is circulant from Lemma 1.1 [1, 
p. 9791]. Let 

J„~ = circ(m 1 ,ni2, . . . , m n ). 



Since the last row of the matrix T n is 



q 1 — Jn 2 J„_2 2J„_ 3 



2J 3 2J 2 



9n 



) * * • ) J 

Qn Qn Qn 



9n - 1 

m 2 = — 

1 f[(l - J„ - .g„)(J„+ 2 - 1) - 4J„J„_ 2 ](1 - J, i+1 )™- 4 

7714 = 



■ (1- J„-.g„)(2" +1 -2)(l 
_2_ Mr-^ , (1 - J n+ i) fc ~ 1 



3 



-1\ (1 — -Wl) 



n— 5 



(2J„)"- 2 
(1 - J„+i) A 



fc=i 



(2J„)» 



"2^ J fc+1 



fc=i 



(1 - J n+l ) k ~ l 

{2J n ) k 



I", 



nil 



-[(1- J„-.g„)(2" +1 -2)(l + (-l)"- 1 ) 



+ 2 J n _2 



2J„ (2J„) 2 J 



+ 2J„_ 3 - 



2t/ n 



Jn+l + (1 - 2J„_i)g„ - 1 



2<?nJ 2 

where g„ = 1 — J„ + 2 2Z^fc=i Jn—k—i ( ) • ^ we rearrange 7715, then 



771,5 = 



1 

9Vi 



(1- J„-.g„)(2" +1 -2)(l + (-l)"- 1 ) (1 



n — 5 



+ 



2 



n+1, 



\n— 4 



<?„ (2J„)™" 3 



(2J„)"- 2 

(1 — Jn+l) n ~ A n T (1 — Jn+l) n 

1 (2J„)"- 3 



(1 - J„+i)™- 4 



(2J„) 



n-2 



(2J„) 



n— 3 



/ l-Z/c+3 - ^fc+2 - 2J fc+ i) 
z — ' > .. ' 



fc=l 



1 (1-J ra+1 )"- 4 

g n (2J„)"" 3 



(1 - Jn+l) 

(2J„) fe 

(l-J„- g „)(2"+ 2 -4) 

(2J„) 2 

2 ( J„_i 



fc-1 



2 ( J n _x + J "- 2(1 ~ J " +l) ) , n is odd 



Then 



1 (1-J n+1 )"- 
3n (2J„)™- i+1 



Jn-2(1-J„ 



(l-J„- g „)(2"+ 2 -4) _ 2 + J„- 2 (l-J„ +1 ) 



n is even. 



(2J„ 



-2 J„ 



■/„-2(l-Jn 



71 is odd 
71 is even. 
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and i = 5, 6, . . . , n. Since the matrix JJ^ 1 is a circulant matrix and its last row 
is known, the proof is completed. H 

A Hankel matrix A = (a^) is an n x n matrix such that ciij = Oj—i j+i. It 
is closely related to the Toeplitz matrix in the sense that a Hankel matrix is an 
upside-down Toeplitz matrix. 

Corollary 6 Let the matrix P n be as in ^). Then 



p-i 



Hi 

H 2 



where H\:=[l, 0, . . . , 0] and H2 is an (n — 1) x n Hankel type in which the first 
row is [J n , J n -i, ■ ■ ■ , Ji] and the last column is [1, 0, . . . , 0] . 

Proof. The matrix P^ 1 is computed easily by applying elementary row 
operations to the augmented matrix [P n : I n \. I 

Theorem 7 The matrix 1 n = circ(jo,ji, ■ ■ ■ ,j n -x) is invertible when n > 3. 

Proof. We show that detQ 3 ) = 104 ^ and detQ 4 ) = -675 ^ by Theorem 
||. Then and are invertible. Let n > 5. The Binet formula for Jacobsthal- 
Lucas numbers yields j n = a™ + /?", where a + /3 = 1 and a/3 = —2. Then we 
have 



r— 1 r— 1 

«(l-q") f3(l-n k k 

= -1 Z7JT + -1 371T> (l-ow ,1-pw 7^0) 



1 — aw" 1 — par 

(a + /3) - (a n+1 + I3" +1 ) + a(3uj k (a n + /3 n ) - 2apw fc 
1 — auj k — (3ui k + al3uj 2k 

l- Jn+1 +2^(2-j n ) 
- ! _ w fe _ 2w2fc . *-M,...,n 

We are going to show that there is no u; fc , fc = 1, 2, . . . , n— 1 such that h(uj k ) = 0. 
If 1 - j n+ i + 2u k {2 - j n ) = for 1 - w fc - 2u 2k ^ 0, then u k = ^+^7 1 ) would 
be a real number. By (||) we would have sin (^f-) = so that uj k = —1 for < 
2^ < 2ir. However u = —1 is not a root of the equation 1 — j'„_|_ 1 +2(2— j n )u = 
(n > 5), a contradiction, i.e., h(uj k ) 7^ for any w fc , where k = 1, 2, . . . , n — 1 
and rt > 5. Thus the proof is completed by [1, Lemma 1.1]. M 

Lemma 8 If the matrix S = (sjj)"J = i is 0/ i/ie form 

2j n -i, i = j 
- 2, j = t + l 
otherwise, 
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then § 1 = {sij)"j2i * s given by 



s / . = J (l+2j„_ 1 )3-+i ' ^ — * 



otherwise. 

,n-2 



Proof. Let B:= S§ 1 = (6^ ) so that 6^ = Ylk=l s ^ s fcj- Clearly 



fe« = (i + 2j„_i) ■ — — i 



1 + 2j n _i 



= 1. 



If j > i, then 



Tl-2 



J T S ii S ij 



fe=l 



(1 + 2j n _ 1 )i- + U + j (1 + 2j n _ 1 )i-+i " 



similar for j < i. Thus SS _1 = / n _2- ■ 
Theorem 9 Let n > 3. TTie inverse of the matrix l n is 



where 
h 



\ x =circ(h 0> h l ,...,h n -x) 



o 
hi 
hi 



J_ / 9j» -18+(10-8j n _ 2 )y n 
2<M V (1 + 2j„_ 1 ) 2 

1 ( 4y n -9 \ 



fa = 



2y„ Vl + 2j n _J 

™- 2 / (2 - j ) fe - 2 

+g(2i fe+1 -i,)^ (1+2 .;_ i)fe _ 

— [((4-j n _ 1 -2y n )(j n+1 -1) 

2y n 

-(2j n _i - i„_ 2 )(l + 2j„_i)) ^ 2 j n \) n - 2 

n—A 



hi := 



J_ / (4 - jn-i - 2y n )(jn+i + 8j» - 2j„_! - 9((-2)» + 1) 
2y„ I (1 + 2j„_ 1 ) 2 

„• (4j»-l ~ 2jn-2)(2 - jn) \ (2 - jj"-'- 1 

■ ?,l+J '^ 1 l + 2i„_! J (l + 2j n _ 1 )™-*+a 

/ory„ = |(4-j, i _i) + iX;fc2 (2jfc - jfc-i) ( nrr^f 1 ) fflrld * = 4,5, ...,n-l. 
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Proof. Let 



1 









2j0-jn-l 

1 



+ kjn-l 



(in— 2 

2j„-i-j„ 



i) + sin-2 



2y„ 





^r(j2 - 2j 3 ) + |i 2 








and G = diag(2, y n ) where y„ = ±(4-j„_i) + ± Efc=2 ( 2 ifc - Jfc-i"* ' 



2-j„ 



i—k 



and 2/„ = ELlJk 



n-1 ^ / l+2jn_ 
2-Jn 



n — fc — 1 

. Then we obtain 
X„3„M„Z„ = G © § 



where Gffi S is the direct sum of the matrices G and S. If T n = M n Z n , then we 
have 

I- 1 = T n (G- 1 ®S- 1 )K n . 

Since the matrix 1 n is circulant, the inverse matrix H^ 1 is circulant from 
Lemma 1.1 [1, p. 9791]. Let 

H" 1 = circ(h , hi,..., h n -i). 

Since the last row of the matrix T n is 



g 1 2 30 - jn- 1 _ 1 2j„_l - j„_2 2j n _2 - jn-3 



2y„ 2y n 2y n 2y Ti 

the last row elements of the matrix l^ 1 are 

1 / 4y„-9 



2y n 



hi = 
h 2 = 



2y n \1 + 2j n _i 
1 



2y n 

n— 2 



(4-jn-i -2y n ) 



(2 -Jn) 



n— 3 



fc=2 



(1 + 2i n _i) 



(2 -Jn) 



n-2 



fc-2 



(l+2j n _!)*-i 
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1 



((4 - jn-i - 2y n )(j n+ i - 1) - (2j„_i - j„- 2 )(l + 2i„_i)) 



(2-J„) r 



(l + 2i„_i)' 



fc=i 



+2^((2 Jfc+1 - Jfc )^ 



(2 - jn) 



fc-1 



(1 + 2j„_ 1 ) fe 



1 f (4 - j n _ i - 2y n )(J n+1 + 8j„ - 2j n _i - 9((-2)" + 1) 

(1 + 2j„_i)a 
(4i„_ a - 2j„_ 2 )(2 - i„) \ (2 - in)""*- 1 



2y« 

-2j„ + j„-i 



1 + 2j„_i 
1 /9j n -18 + (10-8j n _ 2 )</ n 



(1 + 2j„_ 1 )«-*+ 2 



(1 + 2j„-i) 5 



n—k 



and i = 4, 5, . 



where y n = \{A-j n _ x )+\ YZll ^ ~ 3k-i) {^BjT-) 
1. Since the matrix H" 1 is a circulant matrix and its last row is known, the proof 
is completed. ■ 

Corollary 10 Let the matrix K n be as in (Q). Then 



K~ := 



1 
C D 



whe 



C: 



Jn—1 3 71— 2 Jn— 3 



■1± 

2 



(n-l)xl 

and D is the (n— l)x (n— 1) Hankel matrix in which the first row is [J n _i, J n -2, ■ ■ 
and the last column is [J\, 0, . . . , 0] T . 

Proof. The matrix K" 1 is obtained easily by applying elementary row oper- 
ations to the augmented matrix [K n \I n ]. M 



Ji] 



References 

[1] S. Q. Shcn, On the determinants and inverses of circulant matrices with 
Fibonacci and Lucas numbers, Appl. Math. Comput. 217 (2011) 9790- 
9797. 

[2] S. Q. Shcn, On the bounds for the norms of r-circulant matrices with 
Fibonacci and Lucas numbers, Appl. Math. Comput. 216 (2011) 2891- 
2897. 

[3] I. J. Good, On the inversion of circulant matrices, Biometrica, 37 (1950) 
185-186. 

[4] P. J. Davis, Circulant Matrices, Wiley, NewYork, 1979. 



13 



[5] M. Akbulak and D. Bozkurt, On the norms of Toeplitz matrices involving 
Fibonacci and Lucas numbers, Hacettepe J. Math, and Stat. 37 (2) (2008) 
89-95. 

[6] M. Miladnovic and P. Stanimirovic, Singular case of generalized Fibonacci 
and Lucas numbers, J. Koeran Math. Soc. 48 (2011) 33-48. 

[7] G. Y. Lee, J. S. Kim and S. G. Lee, Factorizations and Eigenvalues of 
Fibonacci and symmetric Fibonacci matrices, Fibonacci Quarterly 40 
(2002) 203-211. 

[8] A. F. Horadam, Jacobsthal Representation Numbers, 
http://www.fq.math.ca/Scanned/34-l/horadam2.pdf. 

[9] N. L. Tsitsas, E. G. Alivizatos and G. H. Kalogcropoulos, A recursive 
algorithm for the inversion of matrices with circulant blocks, Applied Math. 
And Comp., 188 (2007) 877-894. 

[10] G. Zhao, The improved nonsingularity on the r-circulant matrices in signal 
processing, Inter. Conf. On Computer Techo. and Development-ICCTD 
2009, Kota Kinabalu, 564-567. 

[11] W. Zhao, The inverse problem of anti-circulant aatrices in signal processing, 
Pacific-Asia Conf. on Knowledge Engineering and Software Engineering- 
KESE 2009, Shenzhen, 47-50. 

[12] D. Serrc, Matrices: Theory and Applications, Springer, New York, 2002. 

[13] R.A. Horn and C.R. Johnson, Matrix Analysis, Cambridge Univ. Press, 
Cambridge, 1985. 

[14] R. Aldrovandi, Special Matrices of Mathematical Physics: Stochastic, Cir- 
culant and Bell Matrices, World Scientific, Singapore, 2001. 



14 



